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We present the two-loop SUSY-QCD corrections to the effective bottom Yukawa couplings within the 
minimal supersymmetric extension of the Standard Model. The effective Yukawa couplings include 
the resummation of the non-decoupling corrections Amj for large values of tg/3. We have derived 
the two-loop SUSY-QCD corrections to the leading SUSY-QCD and top-induced SUSY-electroweak 
contributions to Ami,. The scale dependence of the resummed Yukawa couplings is reduced from 
0(10%) to the per-cent level. These results reduce the theoretical uncertainties of the MSSM Higgs 
branching ratios to the accuracy which can be achieved at a future linear e + e~ collider. 
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The Higgs mechanism [l[ is a cornerstone of the Standard 
Model (SM) and its supersymmetric extensions. The 
masses of the fundamental particles, electroweak gauge 
bosons, leptons, and quarks, are generated by interac- 
tions with Higgs fields. The search for Higgs bosons is 
thus one of the most important endeavors in high-energy 
physics and is being pursued at the upgraded proton- 
antiproton collider Tevatron Q with a centre-of-mass 
(CM) energy of 1.96 TeV, followed in the near future 
by the proton-proton collider LHC @] with 14 TeV CM 
energy. 

The minimal supersymmetric extension of the Stan- 
dard Model (MSSM) requires the introduction of two 
Higgs doublets. After electroweak symmetry breaking 
there are five elementary Higgs particles, two CP-even 
(h, H), one CP-odd (A) and two charged ones (H^). At 
lowest order all couplings and masses of the MSSM Higgs 
sector are fixed by two independent input parameters, 
which are generally chosen as tg/3 = Vi/vx, the ratio of 
the two vacuum expectation values 012, and the pseu- 
doscalar Higgs mass Ma- Including the one- loop and 
dominant two-loop corrections the upper bound on the 
light scalar Higgs mass is M/, < 135 GeV The cou- 
plings of the various Higgs bosons to fermions and gauge 
bosons depend on mixing angles a and /3, which are 
defined by diagonalizing the neutral and charged Higgs 
mass matrices. 

The negative direct searches at LEP2 yield lower 
bounds of M h ,h > 92.8 GeV and M A > 93.4 GeV. The 
range 0.7 < tg/3 < 2.0 in the MSSM is excluded by the 
Higgs searches for a SUSY scale Msusy — 1 TeV at the 
LEP2 experiments 

The dominant genuine SUSY-QCD and SUSY-electro- 
weak corrections to bottom- Yukawa-coupling induced 
processes, as e.g. Higgs boson decays to bb pairs and 
Higgs radiation off bottom quarks, can be derived from 
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where Amb determines the relative corrections to the 
bottom Yukawa couplings. The Higgs couplings are given 
by 
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The leading one-loop corrections Arab to these effective 
couplings can be cast into the form 
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with the scalar function 
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The parameter v = \/v\+v\ = y I/V^Gf is related 
to the Fermi constant Gf and At = v2mj/u2 denotes the 
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top Yukawa coupling. The SUSY-QCD and top-induced 
SUSY-electroweak corrections turn out to be significant 
for large values of tg/3 and moderate or large \i and A t val- 
ues. In order to improve the perturbative result all terms 
of O [{a sl xtgp/M SU SY) n ] and O UXpA t tg/3/Msusr) n ] 
have been resummed in Eq. ^) [7|, |8|. The correction 
Arrib is non- decoupling in the sense that scaling all SUSY 
parameters ,rrig,fi in Eq. |2]) leaves Am& invariant. 
However, its contribution develops decoupling properties 
9\. The corrections Am b contain the strong coupling 
a s (//fj) and the top Yukawa coupling X t (fiR) with signif- 
icant renormalization scale dependences. This leads to 
theoretical uncertainties in e.g. the MSSM Higgs boson 
decay widths and branching ratios of up to 0(10%) @ 
which are larger than the achievable accuracy at a fu- 
ture e + e~ linear collider (ILC) [lcj . In this letter we 
present the two-loop SUSY-QCD corrections to the con- 
tributions Amb of Eq. (UJ) in order to reduce the theoret- 
ical uncertainties to the per-cent level. 
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1: One-loop diagrams of (a) the SUSY-QCD and (b) 
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the top-induced SUSY-electroweak contributions to the bottom 
self-energy with mass insertions corresponding to the correc- 
tions Ami of the bottom Yukawa couplings. The virtual par- 
ticles involve bottom quarks b, sbottom b and stop t squarks, 
gluinos g and charged Higgsinos h . 

Typical SUSY-QCD and SUSY-electroweak diagrams 
that contribute to the bottom self-energy at one-loop or- 
der are displayed in Fig. [1] The leading Amb corrections 
can be obtained from the diagrams in Fig. [1] with off- 
diagonal mass insertions in the virtual sbottom and stop 
propagators in the chiral squark basis 2l|. These mass 
insertions yield a factor Ab/it>2 for the sbottom propa- 
gators and \ t A t V2 in the stop case. One obtains the 
one-loop results of Eq. ||3J) by replacing u 2 
expressing the neutral Higgs component < 
ond Higgs doublet by the mass eigenstates h,H,A [8|. 
These replacements lead to the exact interactions with 
non-propagating Higgs fields, i.e. in the low-energy limit 
of small Higgs momentum [111 ] . This method will be 
applied to the leading two-loop diagrams within SUSY- 
QCD. A typical sample of two-loop diagrams contribut- 
ing to the bottom self-energy is shown in Figs. [5^,b [a 
mass insertion has to be included in all possible ways in 
the sbottom/stop propagators]. 

Dimensional regularization has been adopted for isolat- 
ing the ultraviolet singularities. The bottom momentum 
and its mass have been put to zero while keeping the 
bottom Yukawa coupling Ab finite in the mass insertions. 
All supersymmetric particles as well as the top quark 
have been treated with full mass dependence. The two- 
loop vacuum integrals have been reduced to the two-loop 
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FIG. 2: Typical two-loop diagrams of (a) the SUSY-QCD and 
(b) the top-induced SUSY-electroweak contributions to the bot- 
tom self-energy involving bottom quarks b, sbottom b and stop 
t squarks, gluons g, gluinos g and charged Higgsinos h . 

master integral Tx^{m,\, ms, m^) [121 ] and one-loop one- 
point functions Ao(m) by standard reduction methods 
[131 ] . After adding all two-loop diagrams linear ultravi- 
olet divergences are left over which are absorbed by the 
renormalization of all masses and couplings appearing at 
one-loop order. The heavy masses mt.,m^.,mg appear- 
ing in the propagators have been renormalized on-shell. 
The trilinear coupling A t has been treated in the on-shell 
scheme, too. The strong coupling a s and the top Yukawa 
coupling At have been defined in the MS scheme with 5 
active flavors, i.e. the top quark and the supersymmetric 
particles have been decoupled from the scale dependence 
of the strong coupling a s (fJ,R). Care has to be taken to 
include only the desired order, i.e. 0{a 2 s ^tgf3/MsuSY) 
for Amf D and 0{a s \ 2 t A t tg(3/Ms USY ) for Amf w . In 
this order the trilinear coupling A t only receives a finite 
renormalization at NLO. More details of our calculation 
will be published in [13 ]. 

Since dimensional regularization violates supersymme- 
try by e.g. attributing (n — 2) degrees of freedom to the 
gluons but two degrees of freedom to its supersymmet- 
ric gluino partner, anomalous counter terms have to be 
added in order to restore the supersymmetric relations 
between the corresponding couplings. Since the strong 
coupling factors at one-loop correspond to the Yukawa 
couplings between gluino, sbottom and bottom quark 
states an anomalous counter term has to be introduced 
in order to express all strong couplin g fa ctors in terms of 
the conventional MS QCD coupling [l5j . Moreover, the 
bottom Yukawa coupling Ab in Fig. [TJl determines the 
strength of the Higgs coupling to sbottom states at one- 
loop order which differs from the MS bottom Yukawa 
coupling by a finite amount [l5j |. The total anomalous 
counter term for Amf CD is given by [Ca = 3, Cf = 4/3] 
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For Amf w the situation is similar. The one-loop order 
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involves the bottom and top Yukawa couplings of the 
Higgsino-bottom-stop vertices as well as the top Yukawa 
coupling of the Higgs-stop-stop vertex. Both are shifted 
from the MS couplings by finite amounts 15]. The sum 
of anomalous counter terms for Ami is given by 

The final results have been included in the program HDE- 
CAY 16] which calculates the masses and couplings of 
the MSSM Higgs bosons as well as their decay widths 
and branching ratios. 

The numerical analysis of the corrections Amb and 
their impact on neutral Higgs boson decays is performed 
for the "small a e //" MSSM scenario [ijj as a represen- 
tative case: 

tg/3 = 30, Mq = 800 GeV, M § = 500 GeV (5) 
M 2 = 500 GeV, A b = A t = -1.133 TeV, fi = 2 TeV 

For the Higgs masses and couplings 12211 w e use the RG- 
improved two-loop expressions of Ref. [18|] The bottom 
quark pole mass has been chosen to be Mb = 4.60 GeV, 
which corresponds to a MS mass m^mf,) = 4.26 GeV. 
The strong coupling constant has been normalized to 
a s (M z ) = 0.118. 

The scale dependences of the corrections Am® CD and 
Am* are displayed in Fig. [3] at one- and two-loop or- 
der. The central scale of the SUSY-QCD part Am® CD is 
chosen as the average of the SUSY-particle masses con- 
tributing at one loop, i.e. fio = (mg +mj )2 +mg)/3, and as 
Ho = (mj + mi + /i)/3 for the SUSY-electroweak part 
Amf w . We obtain a significant reduction of the scale 
dependence at two-loop order and thus a large reduction 
of the theoretical uncertainty. Moreover a broad maxi- 
mum/minimum develops at scales of about 1/4 to 1/3 of 
the chosen central scale in contrast to the monotonous 
scale dependences at one-loop order. In the "small a e //" 
scenario the SUSY-QCD corrections are large and posi- 
tive, while the SUSY-electroweak corrections are of mod- 
erate negative size. However, the sign and size of the 
corrections depends on the chosen MSSM scenario. The 
two-loop corrections amount to O(10%) in Am® D and 
a few per cent in Amf w for the central scale choices. 

The branching ratios of the neutral MSSM Higgs 
bosons are depicted in Figs.|4^-c. The bands at one-loop 
order (dashed blue curves) and two-loop order (full red 
curves) are defined by varying the renormalization scale 
between 1 /3 and 3 times the corresponding central scale 
of the SUSY-QCD and SUSY-electroweak parts. We only 
show the two dominant decay modes into bb and t + t~ 
pairs. The uncertainties of the branching ratios reduce 
from O(10%) at one- loop order to the per-cent level at 
two-loop order. The per-cent accuracy now matches the 
expected experimental accuracies at a future linear e + e~ 
collider. 
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FIG. 3: Scale dependence of the corrections Ami, at one- and 
two-loop order: (a) the SUSY-QCD part Amf CD and (b) the 
SUSY-electroweak part Ami in the "small a e ff" scenario. 

Since we have determined the effective resummed 
Yukawa coupling at two-loop order the results will also 
affect all other processes which are significantly induced 
by bottom Yukawa couplings, as e.g. MSSM Higgs radia- 
tion off bottom quarks at e + e~ colliders [lj| and hadron 
colliders [2(J. The two- loop corrections can easily be in- 
cluded in the corresponding numerical programs. 

In summary, the significant scale dependence of 
O(10%) of the NLO predictions for processes involving 
the bottom quark Yukawa couplings of supersymmetric 
Higgs bosons requires the inclusion of NNLO corrections. 
For the corrected Yukawa couplings, we find a reduction 
of the scale dependence to the per-cent level at NNLO. 
The improved NNLO predictions for the bottom Yukawa 
couplings can thus be taken as a base for experimental 
analyses at the Tevatron and the LHC as well as the ILC. 
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FIG. 4: Branching ratios of (a) the light scalar, (b) the heavy 
scalar and (c) the pseudoscalar Higgs boson in the "small 
a e f f " scenario. The dashed blue bands indicate the scale de- 
pendence at one-loop order and the full red bands at two-loop 
order by varying the renormalization scales between 1/3 and 
3 times the central scales given by the corresponding average 
SUSY particle masses. 
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